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ABSTRACT 

Recent observations with the infrared astronomical satehite AKARI have shown that the 
CO2 bands at 4.2 /im in three brown dwarfs are much stronger than expected from the unified 
cloudy model (UCM) based on recent solar C & O abundances. This result has been a puzzle, 
but we now find that this is simply an abundance effect: We show that these strong CO2 bands 
can be explained with the UCMs based on the classical C & O abundances (log^c and logAo), 
which are about 0.2 dex larger compared to the recent values. Since three other brown dwarfs 
could be well interpreted with the recent solar C & O abundances, we require at least two model 
sequences based on the different chemical compositions to interpret all the AKARI spectra. The 
reason for this is that the CO2 band is especially sensitive to C & O abundances, since the CO2 
abundance depends approximately on AqAq — the cube of C & O abundances. For this reason, 
even low resolution spectra of very cool dwarfs, especially of CO2, cannot be understood unless 
a model with proper abundances is applied. For the same reason, CO2 is an excellent indicator 
of C & O abundances, and we can now estimate C & O abundances of brown dwarfs: Three out 
of six brown dwarfs observed with AKARI should have high C & O abundances similar to the 
classical solar values (e.g. log^c = 8.60 and log — 8.92), but the other three may have low C 
& O abundances similar to the recent solar values (e.g. log^c = 8.39 and log^o = 8.69). This 
result implies that three out of six brown dwarfs are highly metal rich relative to the Sun if the 
recent solar C & O abundances are correct. 

Subject headings: brown dwarfs — infrared: stars — stars: abundances — stars: atmospheres — stars: in- 
dividual (2MASS J04151954-0935066, 2MASS J05591914-1404488, 2MASS J15232263+3014562, SDSS 
J053951.99-005902.0, SDSS J083008. 12+482847.4, SDSS J144600.60+002452.0) — stars: low-mass 

1. INTRODUCTION 



'^Also at Department of Astronomy, School of Science, 
The University of Tokyo, Kongo 7-3-1, Bunkyo-ku, Tokyo In our recent work, the CO2 molecule was iden- 

113-0033, Japan tified for the first time in the spectra of brown 

dwarfs observed with the infrared astronomical 



satellite AKARI (jYamamura et al. 1 120101 here- 
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after referred to as Paper I). We tried to inter- 
pret the observed behavior of the CO2 band with 
the use of the model photosphere of brown dwarfs, 
referred to as th e unified cloudy model (UCM; 
TsuiiluOOa I2005I ). In modeling the photospheres 
of brown dwarfs, one problem is how to consider 
the chemical composition, since no direct abun- 
dance analysis is known for brown dwarfs. We 
thought it reasonable to assume a typical compo- 
sition for the disk stellar population such as the 
Sun. However, the solar composition itself experi- 
enced drastic changes in the past decades, and the 
true chemical composition of the Sun is still by no 
means well established. Nevertheless, we thought 
it appropriate to use the latest version of the so- 
lar abundances as possible proxies for the chemical 
abundances in the brown dwarfs. 

In our earlier version of UCMs (lTsuiill20Q2l) . we 
assumed the solar abundances largely based on the 
classical LTE analysis using one dimensional (ID) 
hydrostatic model photospheres and, in particu- 
lar, C & O abundances were log^c = 8.60 and 
IorAq = 8.92 on the sca l e of log Ah = 12 .0 (e.g. 



1997 : lOppenheimer et all Il998l: iGrifhth fc Yelle 



19891 : iGrevesse et all IT99L) . 



Anders fc Grevesse 
At about the same time as we were comput- 
ing our first version of UCMs, a new result for 
the solar C & O abundances based on three di- 
mensional (3D) time-dependent hydrodynamical 
model of the solar photosphere was published 
( AUende Prieto et al.l 2002). Since the classical 
ID model may be too simplified for the real so- 
lar photosphere, this new approach seemed to be 
a useful contribution to the solar abundance anal- 



ysis. The current version of UCMs (|Tsuiill2005[ ) is 
thus based on this new result (\orAc, = 8.39 and 
logAo = 8.69) b y lAllende Prieto et al.l (|2002l) as 



noted elsewhere ( Tsuii et al. 20041) . The new C 



& O abundances are about 0.2 dex smaller as com- 
pared to the classical values referred to above. 

We applied our current version of UCMs to the 
brown dwarfs observed with AKARI in Paper I, 
and we could explain about half of our sample of 
spectra almost perfectly. For the other half of 
our targets, we could explain the overall SEDs 
by this version of UCMs, but we could not ex- 
plain their strong CO2 bands. One explanation is 
that this may be due to an unknown process re- 
lated to CO2, since anomalously strong CO band 
depths have also been explained by a s pecial pro- 
cess now known as vertical mixing (e.g. iNoU et al 



1999[lSaumon et al.ll20o'alLeggett et al.ll2007b[ ) 

However, we happened to try our old version of 
UCMs based on the classical C & O abundances 
and found that the CO2 band appeared to be much 
stronger in the spectra based on the old models 
than on the present models. At the first glance, 
this is rather surprising because C & O abun- 
dances in the old models are only about 0.2 dex 
larger than those in the present models. However, 
we realized immediately that the CO2 abundance 
is extremely sensitive to both C & O abundances 
because the CO2 abundance depends on the cube 
of C & O abundances (^002 °^ ^cAq). We recall 
that the strong dependence of the CO2 abundance 
on metallicity, [Fe/H]0, was previously known by a 
detailed thermochemical analysis of the C, N, and 
O be aring gaseous molecules ([Lodders &: Feglevi 
20021 ). 



The above result demonstrates that at least two 
different series of model photospheres are needed 
for the analysis of the CO2 band observed with 
AKARI. For this purpose, we reconsider our old 
version of UCMs based on the classical C & O 
abundances to represent a case of rather high C 
& O abundances. Our current version of UCMs 
based on the new C & O abundances will serve as 
representing a case of the reduced C & O abun- 
dances compared to the old version of the UCMs. 
An important implication of this result is that the 
metallicity (C & O abundances) in brown dwarfs 
should have a variety of values. 

The interpretation and analysis of the spec- 



historv fe.e. KirkpatricM 2005t Burgasser et al. 


2006a: Leggett et al. 


2007a: ICushing et al. 20081 


Stephens et al.ll2009: 


Yamamura et al. 2010l). and 



the effect of metallicity h as been discussed by some 
authors. For example, Burgasser et al. ( 2006bh 



have measured the strengths of the major H2O 
and CH4 bands in the 1.0-2.5 ^m region in a large 
sample of T dwarfs, and found that the resultant 
spectral indices plotted against spectral type re- 
vealed considerable scatter. Several reasons for 
this result including the effects of dust, gravity, 
and metallicity have been considered, but it ap- 
peared difficult to separate the effect of metallic- 



^The differential iron abundance of a star relative to the Sun 
and defined by [Fe/H]= log(AFc/AH)* - logiApc/Au)Q. 
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ity fro m the remaining paremeters. Leggett et al.l 
lave shown that the effect of metalUcity 
on the SEDs of T dwarfs should be significant, 
but noted that other parameters such as gravity 
can affect the SEDs similarly. This result again 
showed the difficulty in determining metallicity 
uniquely from SEDs. Also, the so-called "blue" L 
dwarfs classified as L subdwarfs have been inter- 
preted to have low i netallicity with [ Fe/HJ from 
-1.5 to -1.0 (e.g. iBurgasser erall l2009l) . but 
those L dwarfs with unusually blue near-infrared 
colors can also be explained by a patchy cl oud 
model ( Folkes et al. I I2OO7I: iMarlev et al.ll2010l) . 

The brief survey outlined above reveals that the 
problem of metallicity in brown dwarfs is still un- 
resolved. In this paper, we will show clear evi- 
dence of metallicity variations in brown dwarfs for 
the first time. In fact, the most important sig- 
nificance of the discovery of CO2 with AKARI is 
that it demonstrated the variations of the C & O 
abundances by at least 50% in brown dwarfs and 
that it provided a means by which to estimate the 
C & O abundances in very cool dwarfs. 

In Paper I, we have analyzed the AKARI spec- 
tra and discussed the basic physical parameters 
of our objects in detail. There we have applied 
the conventional method based on a direct com- 
parison of the observed and predicted spectra. In 
this paper, we examine the results of Paper I by 
a more detailed numerical method in Section 14. 1[ 
and we confirm that the physical parameters de- 
termined in Paper I mostly agree with those based 
on the reduced-chi-square minimization method 
within the estimated errors. A problem, however, 
is that an adequate application of such a rigor- 
ous numerical method requires the input data of 
sufficient accuracy. Unfortunately, the input data 
— our present models of brown dwarfs — are not 
precise enough for this purpose, as discussed in 
Section 4.5 of Paper I. Therefore the numerical 
method does not necessarily provide the best an- 
swer and the traditional fitting method "by eye" 
can still be useful for some cases. For these rea- 
sons, we use the physical parameters determined 
in Paper I and adopt the same approach as Pa- 
per I, eye-fitting, throughout this paper. 




Fig. 1. — The effects of C & O abundances on the 
thermal structure of cool dwarfs. The solid and 
dashed lines are based on the 3D and ID solar 
abundances (see Table [IJ , respectively, a) Toff = 
1500 K, r„ = 1700 K, and \ogg = 4.5. b) T^s = 
1200 K, = 1900 K, and \ogg = 4.5. c) T^s = 
900 K, T„ = Teond, and \ogg = 4.5. 

2. ROLE OF THE CARBON AND OXY- 
GEN ABUNDANCES IN THE UCM 

Our two series of UCMs are referred to as UCM- 
a and UCM-c, and they only differ in C & O abun- 
dances as summarized in Table [TJ The UCM-a se- 
ries is based on the classical C & O abundances, 
which we refer to as ID solar abundances for sim- 
plicity, and the UCM-c series on the new abun- 
dances, which we refer to as 3D solar abundances. 

We first examine the effects of C & O abun- 
dances on the thermal structure of the photo- 
sphere. For this purpose, the models of the UCM- 
c series are taken from our database[l. Since our 
code to compute UCMs has been modified to some 
extent over the last 10 years, we recompute all the 
models of the UCM-a series used in the present 



The numerical details of this ver- 
sion (UCM-c series) are available from 
|http:/ /www. mtk.ioa.s. u-tokyo.ac.jp/~ttsuji/export/ucmLM/"] 
and /ucm/. 
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Table 1: Carbon and Oxygen Abundances in UCMs 



series log Aq " log Aq °- Note on chemical composition 



UCM-a 8.60 8.92 ID solar abundances (e.g. Anders fc Grevesse 1989) '' 

UCM-c 8.39 8.69 3D solar abundances (e.g. AUende Prieto et al. 2002) 



Notes. 

a) The logarithmic abundance on the scale of log Ayj = 12.0. b) Actually, we have applied a slightly 
updated version as summarized in Table 1 of iTsuiil ( 20021) . c) A full listing of the abundances is given in 
[http:/ /www.mtk.ioa.s.u-tokyo.ac.jp/^ttsuji/export/ucm/tables/tablel.dat . 



paper. Therefore the models of the UCM-a and 
UCM-c series are now computed by exactly the 
same code, except for the Rosseland and Planck 
mean opacities, which of course differ according to 
the chemical composition adopted. 

We show a simple comparison of the photo- 
spheric structures of the UCM-a and UCM-c se- 
ries for the cases of T^s = 900, 1200, and 1500 K 
in Figure [H Other parameters such as Tcr and log 
g are chosen to be those actually found for our ob- 
jects (see Table [2|). Inspection of Figure [1] reveals 
that the models of the UCM-a series shown by the 
dashed lines are generally warmer by up to about 
100 K as compared to the models of the UCM-c 
series shown by the solid lines. Since the major 
opacity sources such as CO and H2O are more 
abundant in the UCM-a than in the UCM-c se- 
ries, the blanketing effect due to molecular bands 
should be more effective and hence the models of 
the UCM-a series are warmer than those of the 
UCM-c series. 

Next, we examine the effects of C & O abun- 
dances on the CO2 and other molecular abun- 
dances. We present the abundances of H2O, CO, 
CO2, and CH4 for the case of T^g = 1500 K, 
Tcr = 1700 K, and \ogg = 4.5 for the models of 
the UCM-c ( applied to 2MASS J152322+3014 
in Section and UCM-a (applied to SDSS 
J083008-h4828) series in Figure [2 as the solid and 
dashed lines, respectively. The increased C & 
O abundances result in the increases of CO, CO2, 
and H2O abundances as expected. The increase of 
the CO2 abundance in the UCM-a series is quite 



significant for the reason noted before. On the 
contrary, the CH4 abundance shows a decrease 
in the UCM-a series and this unexpected result 
may be because the direct effect of the increased 
carbon abundance on the CH4 abundance is su- 
perseded by the dissociation of CH4 due to the 
elevated temperatures in the model of the UCM-a 
series (Figure [1]). 

As another example, we show the case of T^s 
= 1200 K, r„ = 1900 K, and \ogg = 4.5 in Fig- 
ure [3l The results are again shown for the UCM-c 
and UCM-a series with the solid and dashed lines, 
respectively. In this case, the effects of the abun- 
dance changes are more pronounced, especially for 
CO2. We will see in Section [3] that the case of the 
UCM-a series is approximately realized in 2MASS 
J055919-1404 in which the CO2 band appears to 
be very strong. 

3. EFFECTS OF THE CARBON AND 
OXYGEN ABUNDANCES ON THE 
SPECTRA OF BROWN DWARFS 

We compare the observed spectra of six brown 
dwarfs with the predicted ones based on the mod- 
els of the UCM-a and UCM-c series in Figure g] 
(a)-(f). The effect of C & O abundances on the 
spectra of brown dwarfs can be seen most clearly 
in a comparison of 2MASS J152322-H3014 and 
SDSS J083008-H4828 shown in Figure Sljc) and 
Figure [4ljd), respectively. These two brown dwarfs 
were found to have nearly the same physical pa- 
rameters (Toff = 1500 K, Tcr = 1700 K, and logy 
= 4.5) but the spectra looked to be quite dif- 
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Fig. 2. — The effects of C & O abundances on the 
molecular abundances in cool dwarfs. The solid 
and dashed lines are based on the 3D and ID solar 
abundances (see Table 1), respectively. The physi- 
cal parameters of the model are: Teg = 1500 K, Ta 
= 1700 K, and log 5 = 4.5. This case may apply 
to 2MASS J152322+3014 (solid line) and SDSS 
J083008+4828 (dashed line) . The units of Pg and 
Pmoi are dyncm^^. The arrow indicates the onset 
of convection. 

ferent (Paper I). In particular, the CO2 band at 
4.2 /Ltm appeared to be much stronger in SDSS 
J083008+4828 than in 2MASS J152322+3014. In 
2MASS J152322+3014, the observed spectrum 
could be accounted for by the model of the UCM- 
c series (Paper I), as confirmed by curve 1 in 
Figure IH^c): Especially, the regions of the II2O 
2.7 /xm and the CO2 4.2 /^m bands as well as the 
Q-branch of CH4 band appeared to be well ex- 
plained by the model of the UCM-c series. On the 
other hand, the predicted spectrum based on the 
model of the UCM-a series shown by curve 2 can 
not explain those features. 



Fig. 3. — The same as Figure[2]but for the physical 
parameters of the model: Tcff — 1200 K, — 
1900 K, and log 5 = 4.5. This case may apply to 
2MASS J055919-1404 (dashed line). 

In SDSS J083008+4828 shown in Figure Hd), 
the observed spectrum could not be accounted for 
by the model of the UCM-c series (Paper I), as 
confirmed by curve 1. On the other hand, the 
strong CO2 band at 4.2 /im can now be explained 
reasonably well by the model of the UCM-a se- 
ries, as shown by curve 2 in Figure |4l[d) . Thus 
the large depression due to the CO2 band turns 
out to be due to the high C & O abundances in 
a LTE model. It is to be noted that this result is 
due to the increase of both C & O abundances. In 
fact, a change of the C abundance alone, for exam- 
ple, produced only minor change on the CO2 band 
strength (Paper I). The large depression over the 
2.7/Ltm region mostly due to II2O can also be bet- 
ter explained with the model of the UCM-a series. 
On the contrary, the Q-branch of the CH4 band 
appears to be weaker with the model of the UCM- 
a than with that of the UCM-c series, consistent 
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Fig. 5. — Observed spectrum of 2MASS 
J041519— 0935 is compared with the predicted 
spectra based on the models of the UCM-a se- 
ries. The best fit is obtained for T^s = 900 K, 
Tci = Tcond, and \ogg — 4.5. a) Effect of Tcff un- 
der the fixed values of Ta = Tcond and log g — 
4.5. b) Effect of log g under the fixed values of 
Teff = 900KandTer = Teo„d. 



with the decrease of the CH4 in the UCM-a com- 
pared to the UCM-c series (Figure [2|). Also, the 
observed CH4 Q-branch indeed agrees better with 
the predicted one based on the model of UCM-a. 
Thus, we conclude that the UCM-a series should 
be applied to SDSS J083008-)-4828 rather than the 
UCM-c series. 

In Figure HJa), we compare the observed spec- 
trum of SDSS J053952-0059 with the predicted 
ones based on the models of the UCM-a and UCM- 
c series. We already know that this spectrum is 
well explained by a model of UCM-c series in Pa- 
per I (curve 1). On the other hand, the observed 
spectrum cannot be explained by a model of the 
UCM-a series of the same parameters (curve 2). 
We obtain more or less similar result for SDSS 
J144600-I-0024, namely the observed spectrum of 
this object can be well explained by the model of 



the UCM-c (curve 1 in Figure HKb)), but not with 
that of the UCM-a series (curve 2). 

In Figure (H^e) , we examine the case of 2MASS 
J055919-1404 in which the observed CO2 band is 
very strong. We could not explain the CO2 and 
CO bands in this object with our UCM-c series 
(Paper I), as confirmed by curve 1. But we can 
now explain the strong CO2 band approximately 
with our model of the UCM-a series, as shown by 
curve 2. This result is fairly consistent with the 
very large increase of the CO2 abundance for this 
model as noted in Figure [3l The fit of curve 2 
based on the UCM-a series can in principle be im- 
proved further by a fine tuning of C & O abun- 
dances. But we defer such a detailed abundance 
analysis to future works and we only note here 
that the strength of the CO2 band is adjustable 
by changing C & O abundances. 

Lastly, we examine the case of 2MASS J041519-0935 
in Figure [3Kf). This is a case in which the fitting 
parameters had to be changed from the result of 
Paper I: We previously found that Tcff = 800 K 
based on the model of the UCM-c series. However, 
the II2O band at 2.7 /zm appears to be too strong 
if we apply the same parameters to the model of 
the UCM-a series. As shown in Figure [5ja), the 
observed spectrum can be explained by the model 
of a higher effective temperature of Tcff — 900 K 
with the model of the UCM-a series. Also, we find 
that the case of logg = 4.5 provides the best fit 
as shown in Figure [HUb) . Thus, we conclude that 
(Teff, T„., \ogg) = (900, T^nd, 4.5) for 2MASS 
J041519-0935 for the models of the UCM-a se- 
ries. We show the best possible predicted spec- 
trum based on the model of the UCM-c series with 
Toff — 800 K and that based on the model of the 
UCM-a series with Teff = 900 K by curve 1 and 
curve 2, respectively, in Figure lU^f). We see that 
the model of the UCM-a series appears to match 
better with the observed spectrum than the model 
of the UCM-c series. 

Finally, we summarize the basic parameters of 
the six brown dwarfs in Table[2l The major change 
from Table 4 of Paper I is to have introduced 
abundance classes in column 3 indicated by the 
UCM series applied : UCM-a and UCM-c means 
that C & O abundances should be close to those 
of the ID and 3D solar abundances, respectively. 
The physical parameters are re-examined with the 
models of the UCM-a series in the same manner 
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Table 2: Basic Parameters from the Model Fittings by Using UCMs 



no. object UCM series'^ Teg (K) r„. (K) \ogg R/Rj- AT^siK)" 



1 SDSS J053952-0059 

2 SDSS J144600+0024 

3 2MASS J152322+3014 

4 SDSS J083008+4828 

5 2MASS J055919-1404 

6 2MASS J041519-0935 



UCM- 


c 


1800 


1800 


UCM- 


c 


1700 


1700 


UCM- 


c 


1500 


1700 


UCM- 


a 


1500 


1700 


UCM- 


a 


1200 


1900 


UCM- 


a 


900 


^cond 



5.5 0.804 -110 

4.5 0.716 -108 

4.5 0.684 -170 

4.5 0.610 -173 

4.5 1.122 -269 

4.5 0.676 -136 



Notes. 

a) The UCM series applied and indicates approximate C & O ab undances (see Tab le [H). b) Radius R relative 
to the Jupiter's radius Rj. c) AT^g = Tcs (empirical values by IVrba et al.l ( 20041 )) — T^s (column 4 in this 
Table). 



as in Paper I. A major change in physical pa- 
rameter is Toff for 2MASS J041519-0935, which 
is changed from 800 K to 900 K as shown in Fig- 
ure [^a). As for other five objects, the physical 
parameters for the models of the UCM-a series 
remain the same as for the UCM-c series. We 
confirm that the overall SEDs based on the mod- 
els of the UCM-a and UCM-c series of the same 
physical parameters agree well (see Figure H]) even 
though some local features due to the CO2 and 
H2O bands differ somewhat. Thus, it is natu- 
ral that the physical parameters mainly derived 
from the fits of the overall SEDs remain the same 
for UCM-a and UCM-c series. The R/Rj values 
for SDSS J083008-I-4828, 2MASS J055919-1404, 
and 2MASS J041519-0935 are changed slightly, 
reflecting the changes of models from the UCM-c 
to the UCM-a series. 

4. DISCUSSION 

4.1. Validity of our "best fit" models 

Fitting of the model spectra to the observed 
spectra is carried out by "eyes" throughout in this 
paper as well as Paper I. The reader might won- 
der whether the eye-fitting can find reliable "best" 
models for various objects. In this subsection we 
assess our eye-fitting by comparing with the nu- 
merical fitting results. 



We evaluate the goodness of fit by the reduced- 
chi-square (hereafter TZ) defined as, 



=1 



/ (A^ - m) , 



(1) 



where fi and Fi are fluxes of the observed and 
model spectra at z-th wavelength grid, respec- 
tively. The uncertainty of the observed flux is in- 
dicated as (7i, and m is the degree of freedom. C 
is the scaling factor that minimizes TZ and is given 

by 



C = 



(2) 



These deflnitions are in principle equivalent to 



"Good ness-of-fit" statistics G by ICushing et al 
'or the equal weight case. 



From our experience in Paper I we know that 
the current UCM cannot fit the observations be- 
yond 4 /zm at least in some objects. Therefore we 
limit the wavelength range for calculating TZ to 
2.64-4.15 /xm (cf. Model spectra are available from 
2.64 /im and CO2 band starts from 4.17/im). 

In Table [3] we list the three best models based 
on the TZ value and the model quoted in Pa- 
per I (by eye-fitting) for each object in our sample. 
The eye-fitting results are indicated in bold-face. 
The models selected by the eye-fitting achieve the 
minimum TZ for SDSS J053951-0059 (L5) and 
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Table 3: The Best Three Models According to the 
Numerical Fitting. 



No. 


Tcff 


logs 


T 

cr 


c 


n 




(K) 




(K) 


(Xl06) 




SDSS J053951-0059 (L5) 


1 


1800 


5.5 


1800 


6.32 


1.127 


2 


1900 


5.5 


1800 


5.82 


1.206 


3 


1900 


5.0 


1800 


6.37 


1.288 


SDSS J144600+0024 (L5) 


1 


2000 


4.5 


1700 


1.53 


0.496 


2 


1900 


4.5 


1700 


1.61 


0.513 


3 


1800 


4.5 


1700 


1.69 


0.553 


10 


1700 


4.5 


1700 


1.79 


0.695 



2MASS J152322+3014 (L8) 



1 


1600 


5.5 


1700 


1.89 


0.675 


2 


1600 


5.0 


1700 


1.91 


0.733 


3 


1700 


5.5 


1800 


1.69 


0.740 


4 


1500 


4.5 


1700 


2.45 


0.793 


SDSS J083008+4828 (L9) 


1 


1600 


4.5 


1800 


3.78 


0.679 


2 


1700 


4.5 


1800 


3.52 


0.711 


3 


1800 


5.0 


1900 


3.44 


0.746 


9 


1500 


4.5 


1700 


4.54 


0.841 


2MASS J055919-1404 (T4.5) 


1 


1200 


4.5 


^cond 


21.8 


0.389 


2 


1200 


4.5 


1900 


20.3 


0.418 


3 


1100 


4.5 


1700 


25.3 


0.482 



2MASS J041519-0935 (T8) 



1 


800 


4.5 


^cond 


29.1 


0.170 


2 


900 


4.5 


^cond 


20.6 


0.173 


3 


900 


5.0 


^cond 


17.2 


0.195 



Notes. 

Models of the UCM-c series are adopted through- 
out as in Paper I. The eye-fitting results quoted 
from Paper I are indicated in bold-face. 



2MASS J041519-0935 (T8; for this particular ob- 
ject we search for the best model among those of 
Tcr = Tcond for the reason outlined in Section 4.3.6 
of Paper I) and the second minimum TZ for 2MASS 
J055919-1404 (T4.5). The difference of 7^ be- 
tween the first and second model for the last case 
is tiny, and the model parameters are within the 
uncertainty we stated in Paper I (±100 K for Tos 
and Tcr, and ±0.5 dex for \ogg). 

For two late-L objects, 2MASS J152322+3014 
(L8) and SDSS J083008+4828 (L9), the differ- 
ences in the model parameters are mostly within 
the uncertainty of cyc-fitting, although the models 
used in Paper I are not included in the numeri- 
cal best three for these dwarfs. Altough \ogg of 
2MASS J152322+3014 differs by 1.0 dex, the eye- 
selected model is in the 4th position in the list, and 
we consider that it is still in the accepted range. 

A significant difference between the two meth- 
ods is found in the L5 dwarf SDSS J144600+0024. 
The numerical fitting suggests Tcff = 2000 K as 
the best, which is 300 K higher than the one se- 
lected by the eye-fitting. The second and third 
are of Tes = 1900 and 1800 K. log 5 and T„ are 
the same in all fours models. In fact we regarded 
such high Tcff values to be unrealistic for an L5 
dwarf, and did not consider t hem in the eye-fitti ng. 
The empirical T^s derived bv lVrba et all (|2004h of 
this object is even low as 1592 K. A key feature 
is the CH4 3.3 /xm band, which appears only in 
the Tcff = 1700 K model. The observed spectrum 
of this source is rather noisy and the detection of 
this band is marginal. If the tiny dip seen near 
3.3 /im in the observed spectrum is actually the 
CH4 band, the eye-fitting results, even if it is not 
perfect, are justified. 

The nature of mid-L to early-T type dwarfs are 
still under debate and their effective temperatures 
might actually spread to higher values. Incom- 
plete atmosphere modeling is another possible rea- 
son. As we discuss in Paper I, current atmosphere 
models for brown dwarfs are still exploratory and 
the UCM is one of such models. There are many 
physical and chemical processes not yet under- 
stood in brown dwarfs. These problems shall be 
attacked and eventually incorporated into future 
model atmospheres, but it is beyond the scope 
of the current paper. In addition, some of our 
AKARI spectra have relatively low S/N. Under 
such circumstances, numerical fitting may not al- 
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ways return a unique and physically reasonable 
solution. On the other hand the eye-fitting would 
give weight to some key features and consider bal- 
ance over the wavelength range. Our goal in this 
paper is to highlight the effects of chemical abun- 
dance in the brown dwarf atmosphere. The com- 
parisons described above well demonstrate that 
the eye-fitting is, even if it is not perfect, useful to 
find reasonable models for our purpose. Therefore, 
we apply the model parameters based on our eye- 
fitting for the six objects including J144600-I-0024 
in the analysis throughout this paper. 

It is noted that the differences between the 
models appear much more prominently over the 
shorter wavelength range especially in J-band, 
even if the spectra in the AKARI wavelength 
range are similar to each other. Consideration of 
near-infrared data such as 2MASS photometry or 
ground-based spectroscopy will enable constrain- 
ing the model parameters better (Sorahana et al. 
in preparation). Such improvements in spectral 
range will also help us to evaluate the goodness of 
the fit in the wavelengths beyond 4 /im. 

4.2. C &: O Abundances in Brown Dwarfs 

The very strong CO2 feature observed with 
AKARI in some brown dwarfs has remained a 
puzzle (Paper I), but we find that this is simply 
due to the effect of C & O abundances. Generally, 
a small change in the chemical composition does 
not have a large effect on the predicted spectra at 
low resolution nor on the thermal structure of the 
photosphere in hotter stars. In fact, this is the 
reason why one dimensional spectral classification 
(e.g. Harvard system) is possible for such stars. 
But in the case of cool stars, a small change of 
the chemical composition is amplified in molecu- 
lar abundances. A drastic example is the spec- 
tral branching of cool giant stars into M, S, and 
C types according to whether the C/0 ratio is 
smaller or larger than unity. In cool dwarfs, the 
change of C & O abundances also produces signif- 
icant effect on the strengths of molecular bands as 
well as on the photospheric structures because of 
the large molecular opacities. 

The results of Section [3] reveal that half of the 
brown dwarf spectra observed with AKARI (i.e. 
SDSS J053952-0059, SDSS J144600+0024, and 
2MASS J152322+3014 ) can be fitted by the pre- 
dicted spectra based on the models of the UCM- 



c series (Paper I). Although the fits are by no 
means perfect, the fits with the predicted spec- 
tra based on the UCM-c series are better than 
those based on the UMC-a series. For this reason, 
C & O abundances in these three brown dwarfs 
should be closer to the recent 3D solar abun- 
dances rather than to the classical ID solar abun- 
dances. On the other hand, the remaining half (i.e. 
SDSS J083008+4828, 2MASS J055919-1404, and 
2MASS J041519-0935) of our sample can be rea- 
sonably accounted for by the models of the UCM- 
a series. Therefore, C & O abundances in these 
three objects should be closer to the classical ID 
solar abundances rather than to the recent 3D so- 
lar abundances. 

Since [Fe/H] of the main sequence stars in the 
Galactic dis k covers the range from —0.8 to -t-0.2 
(e.g. ,Edvardsson et al. I[l993 ). the same metallic- 
ity distribution may apply to brown dwarfs. It is 
certainly only by chance that the brown dwarfs we 
have observed are divided into two groups by C & 
O abundances. Our sample is too small to investi- 
gate the metallicity distribution in brown dwarfs, 
and we hope that this problem can be pursued 
further with a larger sample. 

The problem of the solar C & O a bundances is 



still u nder intensive disc u ssion (e.g. lAvres et a.1 



20061 : ICaffau et al.] I2OO8I : lAsplund et al. I l2009l ) 



Although our problem here is not the solar com- 
position, it is of some interest to know which of 
the proposed solar composition results more real- 
istic for the Sun. If the recent 3D result is more 
realistic for the Sun, three of our sample may have 
about solar composition and the remaining three 
may be about 0.2 dex more metal rich. This means 
that the proportion of metal rich objects with the 
highest [Fe/H] of about +0.2 is quite high in our 
present sample of brown dwarfs. 

4.3. Effects of C &: O abundances on the 
0.9 2.5 /im spectra 

We have shown that C & O abundances have 
significant effects on the 2.5-5.0 /im spectra of 
brown dwarfs, and we now examine their effect 
on the 0.9-2.5 /im spectra. As examples, we com- 
pare the predicted 0.9-2.5 /im spectra of the mod- 
els of UCM-a and UCM-c series for the case of 
r„ = 1700 K, Toff = 1500 K, and \ogg = 4.5, 
together with those for the 2.5-5.0 /im spectra in 
Figure m The major difference between the UCM- 
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a and UCM-c series is that H2O bands at 1.1, 1.4, 
1.9, and 2.7 /im are all stronger in the UCM-a 
(curve 2 in Figure [6]) than in the UCM-c series 
(curve 1), and this is due to a direct effect of the 
increased oxygen abundance (see Figure [2]) ■ 



1) UCM-c/Tc=1700/Te=1500/log g=4,5 

2) UCM-a/Tc=1700/Te=1500/log g = 4.5 








ccy\ A(_^- 






co 









12 3 4 5 

A(;um) 

Fig. 6. — Comparison of the predicted spectra 
based on the models of the UCM-a and UCM-c 
series for the case of = 1700 K, Toff = 1500 K, 
and \ogg = 4.5. The 0.9-2.5 /zm spectrum is con- 
volved with the slit function of FWHM = 500 km 
s~^ (a typical resolution of observed spectra in this 
region), and the 2.5-5.0 /im spectrum with that of 
FWHM = 3000 km s'^ (the resolution of AKARI 
spectra). 

Thus the H2O band strengths depend sensi- 
tively on oxygen abundance, and we may hope to 
determine oxygen abundance from the H2O bands. 
However, we must remember that the H2O band 
strengths also depend on other parameters such 
as Tcr, Teff, and \ogg, and we should encounter 
the same difficulty due to a degeneracy of the pa- 
rameters as noted before by other author s (e.g. 
Burgasser eralll2006bl: lleggett et al.ll2009t ). This 
fact reconfirms the unique role of CO2 as a metal- 
licity (C & O abundances) indicator in brown 
dwarfs. 



5. CONCLUDING REMARKS 

Thanks to the AKARI spectra, we are for the 
first time able to demonstrate that the metallicity, 
more specifically C & O abundances, are impor- 



tant parameters to understand brown dwarf at- 
mospheres. Until now, we have assumed that it 
was sufficient to use one sequence of model pho- 
tospheres based on a representative chemical com- 
position in analyzing low resolution spectra of cool 
dwarfs. We must now admit that such an as- 
sumption is inappropriate, and we should consider 
abundance effects more carefully, especially of C & 
O, in our future analysis of cool dwarfs. Also, we 
cannot use any solar composition for cool dwarfs 
unless this substitution can be justified by a direct 
analysis of the spectra of cool dwarfs. 

It is true that a detailed abundance analysis of 
brown dwarfs is difficult especially with low reso- 
lution spectra, but well defined molecular bands, 
even at low resolution, can be potential abundance 
indicators. We know already that CO2 is a fine 
indicator of C & O abundances. Unfortunately, 
however, CO2 is accessible only from space tele- 
scopes and, moreover, spectroscopic observations 
in the near infrared are mostly neglected by the 
recent space infrared missions. From the view 
point of the study on cool dwarfs (and other cool 
stars), the importance of observing the near in- 
frared spectra (especially between 2.5 and 5.0 /im) 
from space cannot be emphasized too much. 

Although the spectra of brown dwarfs appear 
to be complicated, we are now convinced that the 
spectra of brown dwarfs can basically be under- 
stood on the basis of the LTE model photospheres, 
but only if the chemical composition is properly 
considered. This is a reasonable result for such 
high density photospheres as of brown dwarfs in 
which frequent collisions easily maintain thermal 
equilibrium. Thus the chemical composition is the 
most important ingredient in the interpretation 
and analysis of even low resolution spectra. Now, 
with better observed data for brown dwarfs includ- 
ing those from space, analysis of the spectra and 
abundance determination can be done iteratively 
for brown dwarfs as for ordinary stars. 

Finally, we must remember that a major dif- 
ficulty in the analysis of the spectra of brown 
dwarfs is that we have no model of comparable 
accuracy as for ordinary stars yet. For this rea- 
son, even the accurate numerical method such as 
outlined in Section 14.11 cannot be infallible. In 
fact, we have no model reproducing all the ob- 
servable features correctly, and the model found 
by the numerical method as well as by the eye- 
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fitting method may prove ineorreet even if they 
are relatively satisfactory among the models cur- 
rently available. Within this limitation, we hope 
that our main results on the differential effects of 
C & O abundances are relatively free of present 
brown dwarf model uncertainties. 

We thank an anonymous referee for critical 
reading of the text and for invaluable suggestions 
regarding the method of analysis of the spectra 
of brown dwarfs. We are grateful to Dr. Poshak 
Gandhi for his careful checking of the manuscript 
and many suggestions to improve the text. This 
research is based on observations with AKARI, a 
JAXA project with the participation of ESA. We 
acknowledge JSPS/KAKENHI(C) No.22540260 
(PI: I. Yamamura). 
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(a) SDSS J053952-0059 (L5) 

1 ) UCM-c/Tc=1 800/Te=1 800/logg = 5.5 

2) UCM-a/Tc=1800/Te=1 800/logg = 5.5 
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(e) 2MASS J055919-1404 (T4.5) 
1 ) UCM-c/Tc=1 900/Te=1 200/logg=4.5 
- 2) UCM-a/Tc=1900/Te = 1200/logg = 4.5 
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(b) SDSS Jl 44600 + 0024 (L5) 
1 ) UCM-c/Tc=1 700/Te= 1 700/logg=4.5 
- 2) UCM-a/Tc=1700/Te=1700/logg = 4.5 




-(d) SDSS J083008 + 4828 (L9) 

1 ) UCM-c/Tc=1 700/Te= 1 500/logg = 4.5 

2) UCM-a/Tc = 1700/Te=1500/logg = 4.5 
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(f) 2MASS J041 51 9-0935 (T8) 

1 ) UCM-c/Tc=Tcond/Te=800/logg=4.5 

2) UCM-a/Tc=Tcond/Te=900/logg = 4.5 
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Fig. 4. — Spectra observed with AKARI are compared with the predicted spectra based on the models of 
the UCM-c (curve 1, blue) and UCM-a (curve 2, green) series. The UCM-a series is based on the classical 
solar C & O abundances (ID abundances, see Table [IJ and the UCM-c series on the more recent solar C 
& O abundances (3D abundances, see Tabled]). Note that the brown dwarfs shown in (a), (b), and (c) are 
relatively well fitted by the models of the UCM-c series while those in (d), (e), and (f) are better described 
by the models of the UCM-a series. 12 



